Fourth-order magnetic anisotropy and tunnel splittings in Mni2 from 

spin-orbit-vibron interactions 
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From density-functionabtheory (DFT) based methods we calculate the vibrational spectrum of the 
Mni20i2(COOH)i6(H20)4 molecular magnet. Calculated infrared intensities are in accord with 
experimental studies. There have been no ab initio attempts at determining which interactions 
account for the fourth-order anisotropy. We show that vibrationally induced distortions of the 
molecule contribute to the fourth-order anisotropy Hamiltonian and that the magnitude and sign of 
the effect (-6.2 K) is in good agreement with all experiments. Vibrationally induced tunnel splittings 
in isotopically pure and natural samples are predicted. 
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The possibility that vibrationally induced modifica- 
tions of the spin-orbit interaction affects magnetic reori- 
entation barriers has not been previously investigated. 
A simple model captures the physics but quantitative 
determination of the coupling constants requires compu- 
tationally demanding DFT methods. A model Hamil- 
tonian for a single uniaxial anisotropic spin coupled to 
a one-dimensional harmonic oscillator is given by iJ = 
7,,52 + i(p2 + ^2g2) +QJ2abl'abSaSb whcrc a and b 
are summed over x, y, and z, and 7^^ — ^^r)- Complet- 
ing the square shows the diagonal energy of a harmonic 
oscillator {\4>)) and spin {\SM)) product state as a func- 
tion of the magnetic quantum number (M) is given by: 



(1) 



with A = S{S + 1)(7^^ + lyy)/2 and B = -f^^ - (7^^ + 
^yy)/2. For this simple case, the interaction between 
vibrational and spin degrees of freedom always acts to 
further stabilize each M state but the energy splitting 
between the \M\ = S and M — states may be ei- 
ther enhanced or reduced depending on the derivatives 
of the 7 matrix. A more detailed analysis of this Hamil- 
tonian, given below, shows that it can also connect M 
and M ± 4 levels. Further, in a real system, the fre- 
quency, w, depends on nuclear masses so isotope effects 
can lead to small changes in the barrier or, as discussed 
below, larger tunnel splittings in the anisotropy Hamil- 
tonian. Both the intrinsic 4th-ordcr and isotope effects 
can mediate resonant tunneling of magnetization which 
is of great interest from the standpoint of quantum me- 
chanics at the mesoscale. In addition, the magni- 
tudes of terms in the anisotropy Hamiltonian determine 
the suitability of nanoscale particles for use in magnetic 
memory and quantum-computing applications. [^,^ In 
general, purely electronic contributions to higher-order 
magnetic anisotropy scale as [l/(2c^)]'' (c = speed of 
light) but the spin-orbit-vibron contribution scales as 
[l/uj"^] X [l/(2c^)]''. The large prefactor suggests 



this term could be dominant in high-symmetry bulk sys- 
tems where second-order effects vanish by symmetry. 

Recent experiments on Mni20i2(COOR)i6(H20)2 
molecular nanomagnets, commonly referred to as Mni2- 
Ac, [Q-^jjoj provide an ideal system for understanding how 
vibrational degrees of freedom may enhance the mag- 
netic barriers. Many researchers have noted that the tun- 
neling dynamics may be mediated by conventional spin- 
phonon interactions, pof Further, recent experiments by 
Sushkov et al. show some strong variation in infrared 
(IR) spectra as a function of applied field suggesting 
a coupling between the spin and vibrational degrees of 
freedom. Accurate experiments have determined that 
this easy-axis uniaxial molecular magnet has a second- 
order anisotropy parameter oi D — —0.56 K and that 
the fourth-order contributions to the anisotropy Hamil- 
tonian increase the barrier from 56 K to approximately 
65 K. 1^,^,^ Pederson and Khanna have combined DFT, 
spin-orbit coupling and second-order perturbation the- 
ory to calculate the second-order parameter (D) and find 
D=-0.557 K in excellent agreement with experimental 
measurements. [0 An estimate of the entirely electronic 
spin-orbit-induced fourth-order anisotropy may be ob- 
tained by comparing the barrier computed with spin- 
orbit treated within second-order perturbation theory to 
that computed using exact diagonalization. Within this 
approach we find a fourth-order contribution to the bar- 
rier on the order of 1 K and that it acts to reduce rather 
than to enhance the barrier. 

To determine if vibrational coupling is an important 
contribution to magnetic anisotropies in Mni2-Ac we 
have performed accurate DFT based calculations on 
a single molecular unit. The Perdew-Burke-Ernzerhof 
generalized-gradient approximation has been used for 
all calculations. ||l2| Using the methods discussed in 
Ref. p^ , the vibrational frequencies, vibrational eigen- 
vectors, and IR and Raman spectra have been calculated. 
Because this requires a large number (163) of unsym- 
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metrized molecular calculations we have used a Troullier- 
Martins type [ p^|JT5[ | pseudopotential for the O and C 
atoms and have treated the H and Mn atoms within an 
all-electron method. Gaussian basis sets have been fully 
optimized for each atom using the methods of Ref. . 
To determine the anisotropy Hamiltonian and the deriva- 
tives with respect to each normal mode, we calculate the 
anisotropy Hamiltonian H = J^ab^yab^aSb for each of 
the 163 inequivalent vibrational displacements using the 
method of Ref. The coefficients jxy for an arbitrary 
configuration of atoms are determined from matrix ele- 
ments of the spin-orbit-coupling operator sandwiched be- 
tween all pairs of occupied and unoccupied Kohn-Sham 
wavefunctions (squared) and appropriate energy differ- 
ence denominators. The wavefunctions and thus the 
anisotropy Hamiltonian depend on the geometry of the 
molecule. Derivatives of the anisotropy matrix with re- 
spect to the j"* normal mode (i.e. may be deter- 
mined using a finite-difference approach. We have ascer- 
tained that the partial-pseudopotential-based anisotropy 
Hamiltonian for the equilibrium geometry reproduces the 
all-electron anisotropy Hamiltonian of Ref. [0 . To test 
the accuracy of the vibrational frequencies we compare 
our calculated IR spectra, plotted in Fig. 1, directly to 
experimental measurements of Sushkov et al. 

The experiments measure the IR absorption of Mni2 
crystals suspended in paraffin pellets at wavenumbers 
ranging from 30 cm^^ to 70 cm~^ and from 140 cm~^ 
to 650 cm^^. The experimental absorption peak at 
38 cm~^ has a clear corresponding feature in our cal- 
culations at 63 cm~^. The structure in the 140 cm^^ 
to 300 cm~^ range is well reproduced by the calcu- 
lations: The small experimental peaks at 150 cm~^, 
170 cm~^, and 200 cm"'^ correspond to simulated fea- 
tures at 144 cm~^, 170 cm~^, and 201 cm~^, respec- 
tively, and the relative intensities are moderately well re- 
produced. The small peak in the simulation at 180 cm~^ 
can be tentatively identified in the experiment as a small 
peak at 185 cm~^. The structure between 215 cm~^ and 
235 cm""'^ in the experiment, consisting of a large peak, 
a small peak, and a large peak with a shoulder at high 
frequency, has a clear analogue in the simulation results 
between 230 cm~^ and 260 cm^^ with similar relative in- 
tensities. The intense peak at 255 cm^^ corresponds to 
the simulation peak at 275 cm~^, with the smaller exper- 
imental peak at 270 cm~^ appearing in the simulation as 
a shoulder at 260 cm~^. The double peaked structure ob- 
served in experiment at 284 cm^^, which has significant 
activity in a magnetic field, corresponds to the simulated 
peaks at 302 cm~^, 313 cm~^ and 316 cm~^. The fol- 
lowing experimental triplet, consisting of a large peak 
at 300 cm~^ followed by two small peaks at 320 cm~^ 
and 340 cm~^, is present in the simulation, although at 
the resolution of the graph in Fig. 1 the two peaks at 
326 cm~^ and 329 cm~^ overlap, and the third peak ap- 



pears at 343 cm~^. The small peak at 360 cm~^ and in- 
tense peak at 375 cm"-'^ are reproduced in the simulation 
with opposite relative intensities (or interchanged fre- 
quency order) in the double peaked structure at 360 cm~^ 
to 368 cm"-'^. The intense peak at 410 cm~^ bracketed by 
two smaller peaks at 395 cm^^ and 415 cm^^ appears in 
the simulation between 412 cm~^ and 439 cm^^. The two 
small peaks in the calculated IR spectrum at 385 cm~^ 
and 395 cm~^ are not clearly visible in the experimen- 
tal data, although they could correspond to the small 
feature at 382 cm~^ to 390 cm~^ in the gap between 
the prominent double peak and triple peak structures. 
In the highest frequency range measured by experiment 
agreement with the simulation is still good, although the 
relative intensities of the peaks are less accurately repro- 
duced. The two faint experimental peaks at 465 cm~^ 
and 495 cm~^ probably correspond to the two relatively 
prominent peaks at 460 cm~^ and 476 cm^^. The five 
peaks between 510 cm~^ and 570 cm"-'^ correspond to the 
five peaks between 500 cm~^ and 545 cm"-*^. Although 
the intensities relative to the previous two peaks are lower 
in the simulation, opposite to experiment, the relative 
intensities within the five peak structure are good. In 
both the experiment and simulation the second and fifth 
peaks are most prominent. The intense double peak at 
605 cm~^ to 610 cm~^ clearly corresponds to the sim- 
ulated double peak at 564 cm~^ to 571 cm^^. There is 
no comparably intense feature in the simulation that cor- 
responds to the highest frequency experimental peak at 
640 cm^^, although several less intense peaks are present 
in the correct frequency range. 

The decomposition into contributions from different 
structural elements in the Mni2 molecule shows some 
general trends as a function of frequency, but no sharp 
variation from mode to mode. Below 500 cm~^, the 
main contributions are from the Mn and formic acid. 
The water molecules contribute mainly for modes below 
about 350 cm^^, while the anionic oxygen atoms con- 
tribute mainly above this frequency. Above 500 cm^^ the 
Mn contribution grows relative to the formic acid contri- 
bution, with particularly large Mn and anionic oxygen 
weights in the 564 cm~^ to 571 cm^^ double peak. 

We now discuss our calculations on the vibrational 
contributions to fourth-order anisotropy. The coupled 
spin/vibron Hamiltonian is given by: 

H^Y^^abSaSb (2) 

ah 




The above expression may be analyzed either classically 
or quantum-mechanically to determine how vibrational 
coupling affects the magnetization barrier. As demon- 
strated below the magnetization barrier is enhanced by 
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approximately 5.6 K classically and 6.2 K quantum me- 
chanically since the quantum mechanical expection value 
of is 10 % larger than the classical value for S=10. 
However, to address questions related to isotopic effects 
and changes in tunnel barriers a quantum mechanical 
analysis is required. 

The simplest level of approximation is to take a prod- 
uct wavefunction of the form |^) — Uj \SM), with 
\(f)j) a harmonic oscillator wavefunction and \SM) a spin- 
wavefunction, and determine the diagonal energy of this 
state. The resulting energy E{M) = {SM\H\SM)] is 
given by: 



i?(Af)=^7a6C/'^''(M)+Ej' 



(3) 
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{SM\SaSb\SM)] The 4th-order barrier enhancement 
may then be estimated by evaluating the second line of 
Eq. (3) for the M = ±10 and M = states and taking 
the difference. This energy difference is 6.2 K. 

To compare directly to the experimental parametriza- 
tion of Barra et al. which uses 4*''-order Stevens op- 
erators, it is instructive to expand both our calculated 
and their fitted 4th-order expression in terms of orthog- 
onal cubic polynomials of degree 4. Because of the 
S4 symmetry of the molecule, it is generally possible 
to write the fourth-order energy according to: = 
Al+{Ai{^)[S^(Sl-S^/?,)]+A2{A)[?,S^+?,bSt-iOS^Sl] + 
B^{A)[St + 54 - QSlSl] + B^mS^SyiSl - Sl)])/IQ\ In 
Table I, we compare the values of the expansion coeffi- 
cients for each representation. The O4 Stevens operator 
can be expanded into Legendre polynomials of degree 2 
and 4 respectively. 

An interesting feature, present in both the experimen- 
tal and theoretical expansions, is that for the most part, 
the 4*'' order terms exhibits nearly the same angular 
variation as the second-order Hamiltonian. For exam- 
ple the entirely diagonal 4"* order S'^{Sl — 5*^/3) term 
accounts for almost all of the barrier enhancement. This 
feature is very important from the standpoint of sharp 
tunneling transitions. It ensures that magnetic-field in- 
duced alignment condition aligns multiple pairs of states 
simultaneously which opens multiple tunneling paths. 
The coefficient ^2(4) changes the diagonal energies by 
^2(4) X Af^/lO"* and is partially responsible for a broad- 
ening of the magnetic-field alignment condition but does 
not cause any tunnel splittings at zero-field. The Si (4) 
coefficient corresponds to the Bf Stephenson coefficient, 
and leads to observable tunnel splittings in the M = ±2 
manifold and smaller tunnel splittings in the M — ±4 
manifold. Our isotopically pure calculations show tun- 
nel splittings of the M = ±2, M = ±4, and Af± = 6 
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FIG. 1. Calculated total and IR vibrational density of 
states for the Mni2 molecule. In addition we have projected 
the IR active density of states onto Mn, 0^~, COOH, and 
H2O to show the origin of the IR spectrum. 



manifolds of approximately 10^^, 7 x 10^^ and 6 x lO^^*' 
K respectively but do not split the odd M states since 
coupling is between M and M±4. By randomly changing 
a single atomic mass by 1 amu, we find that the odd-M 
states are split with a tunnel splitting of the M — ±1 
states on the order of 10""* K. The isotopic effects exper- 
imentally observed by Wernsdorfer et al. ||l^ in Fes may 
be partially due to the mass effects identified here. 

In Fig. 2 we show the barrier as a function of the num- 
ber of vibrational modes that are included. When none of 
the vibrational modes are included the barrier reproduces 
the earlier calculation of Ref . ||l^ . As expected qualita- 
tively, most (85 %) of the fourth-order barrier is asso- 
ciated with the frequencies in the 100-500 cm"-'^ range 
that correlates with Mn vibrations. The upper panel of 
the figure shows that there are a 5-10 modes which ac- 
count for about 80-90 % of the fourth-order anisotropy 
barrier. For example, single modes at 189, 270, 280, 
and 465 and 1496 cm"! contribute 0.8 K, 1.6 K, 0.8 K, 
0.25 K and 0.45 K respectively to the 4th-order barrier. 
These modes and all other modes that contribute visible 
spikes in Fig. 2 share one simple trait. They are optically 
silent to IR absorption and transmission but show strong 
Raman intensity. The calculated Raman intensities for 
these modes are shown in Table II. 

Recently, Sushkov et al. Q have performed experi- 
ments on Mni2-Ac and have shown that the IR trans- 
mittances at 284, 306 and 409 cm^^ exhibit strong de- 
pendencies when magnetic fields are applied to the crys- 
tal and that the dependence is particularly strong for 
the mode at 284 cni^^. This suggests that vibrations in 
this energy range are associated with the spin carrying 
Mn ions as confirmed from the projected DOS in Fig. 1. 
While the Raman modes are responsible for the forma- 
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FIG. 2. Evolution of 4th order anisotropy barrier (A(4)) as 
a function of the number of modes included. Also shown is 
the total and Mn-projected calculated IR and Raman spectra. 
Large jumps in the barrier are due to strong Raman modes, 
as discussed in the text. 

tion of the 4th-order barrier, the strong contributions due 
to Mn motion at 270 and 280 cm^^ is quahtatively con- 
sistent with the strong field-dependent IR dependencies 
observed experimentaUy in this energy range. 

TABLE L Calculated Raman intensities for modes that 
are primarily responsible for the formation of the vibrational 
4th-order magnetic anisotropy. Units of Raman intensities 
are respectively A^/amu. 



llI (cm ^) Raman 4th-ordcr shift (K) 

189 297 0.8 

270 508 1.6 

281 24 0.8 

465 263 0.25 

1496 5400 0.45 



TABLE IL Fourth-order anisotropy Hamiltonian as deter- 
mined from experiment, DFT plus vibration-spin coupling, 
and DFT without this coupling. For simplicity we have reex- 
panded the representation of Barra et al. in terms of orthog- 
onal cubic polynomials rather than the Stephens polynomials 
used in their work. All numbers are given in KxlO*. 





Ai(4) 


A2(4) 


Bi(4) 


B2(4) 


Experiment 


-8.35 


-0.334 


-0.43 


0.000 


Vibrational 


-5.58 


-0.008 


-0.01 


-0.015 


Electronic 


0.68 


0.0005 


-0.002 


0.004 



While the calculated 4 -order contribution to the tun- 
neling barrier of 6.2 K is in close agreement with the ex- 
perimental values of 7-10 K, it is indeed a small number 
so it is appropriate to consider other vibrational effects 
that might be nonnegligible. 

To determine possible effects of methyl termination of 
the carboxyl groups, we have changed the mass of the 



H on the formate groups from 1 to 15 (the mass of a 
CH3). This further increases the 4th-order anisotropy 
from 6.2 K to 7.3 K and leads to even better agree- 
ment with experiment. There could also be terms due to 

(P^/dQ"^. However, such interactions effectively change 

2 

the vibrational spring constant and scale as rather 

4 

than -^^-^ where /i is the total magnetic moment. Be- 
cause of the large masses involved (even for hydrogen) an 
order of magnitude estimate suggests that the total bar- 
rier would change by less than 0.01 K from this type of ef- 
fect. Further, in contrast to terms discussed above, these 
contributions add as amplitudes and could partially can- 
cel one another. The spring-constant terms directly mod- 
ify the 2nd-order 7 matrix allowing for isotopic induced 
symmetry breaking. While potentially unimportant from 
the standpoint of barrier formation, this effect leads to a 
type of locally varying second-order transverse anisotropy 
required to explain experiments (See Ref. ||T^.) We are 
in the process of studying this and other possibilities. 

To summarize, we have performed accurate DFT- 
based calculations on the Mni2-Ac molecule to deter- 
mine whether vibron-spin coupling could be responsible 
for part of the 4th-order anisotropy Hamiltonian. Our re- 
sults suggest that vibron-spin coupling accounts for most 
of the effect and that it is significantly larger and a differ- 
ent sign than the 0(l/c'*) terms that arise from an exact 
diagonalization treatment of spin-orbit coupling. 

This work was supported in part by ONR grant 
N0001400WX20111 and N0001401WX31303. 
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